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NOTES ON THE TWO-DIMENSIONAL FRACTIONAL 
BROWNIAN MOTION 

By Fabrice Baudoin and David Nualart 1 

Toulouse University of Technology and Universitat de Barcelona 

We study the two-dimensional fractional Brownian motion with 
Hurst parameter H > |. In particular, we show, using stochastic cal- 
culus, that this process admits a skew-product decomposition and 
deduce from this representation some asymptotic properties of the 
motion. 



1. Introduction. We consider a complex fractional Brownian motion started 
at z = x 1 + ix 2 G C*, 

B t = B\ + iB 2 , t>0, 

where B 1 and B 2 are two independent fractional Brownian motions with 
the same parameter H G (0, 1). That is, for i = 1,2, B l = {B\,t > 0} is a 
Gaussian process with mean E(i?^) = x l and covariance 

D U \ 1 I 2H i ,2H i, |2ff\ 
K{t,S) = 2{S +t —\t — s\ j. 

Notice that, in the case H = ^, B is a planar Brownian motion. In what 
follows, we shall always assume that H > ^. In this case, it is known that B 
is a transient process (see [20]). 

The aim of this paper is to use the stochastic calculus for the fractional 
Brownian motion which has been recently developed by several authors (see, 
e.g., [1, 2, 4, 11]) in order to obtain geometric properties of this motion as 
it has been made in the case of the planar Brownian motion (see, e.g., 
[6, 8, 12, 13, 17, 18]). The key point of our study is an analogue of the 
celebrated skew-product decomposition of the two-dimensional Brownian 
motion. Precisely, we show that we can write 

/ f* dB s \ 
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where the integral of the right-hand side can be understood not only path- 
wise, but also as a local divergence. We shall see that it is then possible to 
apply the ergodic theorem to study the asymptotics of the integral 



(1.1) 



dB 



B„ 



whose real part (resp. the imaginary part) explains the radial part of B 
(resp. the angular part of B). 

In that respect, the study of the two-dimensional fractional Brownian 
motion with Hurst parameter H > ^ could seem simpler than the study of 
the planar Brownian motion, for which it is not possible to apply directly 
the ergodic theorem. Nevertheless, for the fractional Brownian motion, we 
shall see that the study of the windings is much more difficult because the 
integral (1.1) is not a time-changed fractional Brownian motion. 

2. Ito's formula for holomorphic functions. Before we study the skew- 
product decomposition, we need to develop stochastic calculus for the two- 
dimensional fractional Brownian motion. In particular, we shall need a two- 
dimensional Ito's formula. 

Let Ti be the Hilbert space defined as the closure of the set of step func- 
tions on R + with respect to the scalar product 

( 1 [o,t]i 1 [o,s])w = R(t,s). 

For i = 1,2, the mapping 1 r ^ — ► B\ — x 1 can be extended to an isome- 
try between Ti and the first chaos H\ associated with B l . We denote this 
isometry by (p If v? G W® 2 , we set B(tp) = (B 1 ^ 1 ), B 2 (ip 2 )). 

Let S be the set of smooth and cylindrical random variables of the form 

(2.1) F = f(B(<p 1 ),...,B(<p n )), 

where n > 1, / G C£°(]R 2n ) (/ and all its partial derivatives are bounded), 
and tfi G TL® 2 . The derivative operator D of a smooth and cylindrical random 
variable F of the form (2.1) is defined as the 7^® 2 -valued random variable 



where Z = (B(ipi), . . . , B(p n )), and we make use of the notation f(z±, . . . , z n ) 
and Zj = (xj,yj), for j = 1, . .. ,n. The derivative operator D is then a clos- 
able operator from L p (0) into L P {Q,;TL® 2 ) for any p > 1. For any p > 1, the 
Sobolev space B 1 ^ is the closure of S with respect to the norm 



\l P = E\F\ p +E\\DF\\* i9 



2 ■ 
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In a similar way, given a Hilbert space V, we denote by D 1 ' p (V r ) the corre- 
sponding Sobolev space of ^-valued random variables. 

The divergence operator 5 is the adjoint of the derivative operator, defined 
by means of the duality relationship 

(2.2) E(FS(u)) = E(DF, u) H ®2 , 

where u is a random variable in L 2 (Q;TC® 2 ). That is, we say that u belongs 
to the domain of the operator 5, denoted by DomJ, if the mapping F t— > 
E(DF,u)f^S)2 is continuous in L 2 (Q). If u £ Dom<5 is a two-dimensional 
stochastic process, we will make use of the notation 



/•OO /"OO 

*(«)=/ u l s 5B]+ u 2 JB 2 . 
Jo JO 



The space D 1,2 (H® 2 ) is included in Dom5. For an element it in B 1 ' 2 (7i® 2 ), 
we have 

(2.3) E{5{u)f = E\\u\\ 2 Hm + E{Du,{Du)*) Hm . 

We recall that "H contains the space of real-valued functions (p on M + such 
that 

rOO /"OO 

|^(t)||(^(s)||t - s\ 2H ~ 2 dtds < oo. 



10 Jo 

The following result has been proved in [2]. 

Proposition 1. Let u = (ut)t>o be a two-dimensional stochastic process 
in the space H 1,2 {TL® 2 ) such that, for some T > 0, 

E( f [ \u t \\u s \\t -s\ 2H ~ 2 dt ds ) < oo, 



o Jo 



E / \D t u e \\D sUri \\t - s\ 2H - 2 \6 - rj\ 2H - 2 dtdudO dr] < oo 

J\0.T} 4 



and 

rT i-T 



[ I \D s u t \\t - s\ 2H 2 dsdt<oo, 
Jo Jo 



io Jo 

a.s. Then the symmetric integral defined as the limit in probability 
/ ut dB t = lim(2e)~ 1 / u s (B, s+e)AT - B (s _ £)v0 ) ds 

JO elO Jo 

exists and we have 

(2.4) / u t dB t = 5(u)+a H f I {D l s u\ + D 2 u 2 )\t - s\ 2H ~ 2 ds dt, 
Jo Jo Jo 

where an = H(2H — 1). 
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We denote by dJ^W® 2 ) the set of ?^® 2 -valued random variables u such 
that there exists a sequence {(n n ,u n )},n> 1} C T x D 1 - 2 ^® 2 ) such that 
Q n ]Q a.s. and u = u n , a.e. on [0,T] x O n . We then say that {(ft n ,u n )} 
localizes u in D 1 - 2 (W® 2 ). If u £ ©^(W® 2 ), by the local property of the 
divergence operator, we can define without ambiguity 5(u) by setting 

5{u)\nn = d(u n )\un 

for each n > 1, where {(f2 n ,u n )} is a localizing sequence for u in D 1 ' 2 ('H® 2 ). 

If / is a real-valued function on M. 2 twice continuously differentiable, then 
for each t > 0, V/(-B.)l( ,t) belongs to D 1 1 C J 2 ('H® 2 ) and the following version 
of Ito's formula holds (see [2]): 



On the other hand, we will denote by /q u s dB s the Riemann-Stieltjes 
integral of a stochastic process u% whose trajectories are 7-Holder continuous, 
provided 7 > 1 — H (see [21]). 

If / : C — > C is a twice continuously differentiable function in the variables 
(x,y) — ► f(x + iy), formula (2.5) leads to 



f(B t ) = f(z ) + J q J-{B s )5B] 



(2.6) 

+ f%{B s )5B 2 s + H l\f{B s )s 2H - 1 ds, 
Jo oz Jo 

Proposition 2. Let f:C — >C be a holomorphic function, then 

(2.7) f(B t )=f(z )+ff'{B s )5B s 

Jo 

(2.8) = f(z )+ [ f'(B s )dB s , 

Jo 

where f'(z) = j£. 

Proof. Equation (2.7) follows from (2.5) and the fact that §= = and 
A/ = 0. Equation (2.8) follows from the change of variables formula for 
integrals with respect to Holder-continuous functions. □ 
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Remark 3. Unlike two-dimensional Brownian motion, two-dimensional 
fractional Brownian motion is not conformal invariant. It means that, in 
general, /q f'{B s ) 5B S is not a time-changed fractional Brownian motion. 
Nevertheless, let us recall that, in [3], it is shown that, for a determinis- 
tic function, g which satisfies some regularity assumptions and for a one- 
dimensional fractional Brownian motion B H with Hurst parameter H > \ , 
we can write a decomposition 

f\(s)dBf = p^ g{s)1/Hds + A t , t>0, 

where At is a bounded variation process and H a fractional Brownian 
motion with Hurst parameter H. Hence, a natural question is to ask if 
such a decomposition exists (with the suitable time-change) for the two- 
dimensional process f(Bt) or not, but so far there is no answer yet to this 
question. 

As an immediate corollary of the previous proposition, we also deduce the 
following: 

Corollary 4. Let a = f(z)dz be a holomorphic one-form on C, then 




where Jg(ot) a denotes the integral of a along the trajectories of B. 

Remark 5. We recall that we can define the integral of a holomorphic 
one-form on any continuous path. 

Remark 6. The previous equality is interesting because it relates ob- 
jects of different natures. 

From now on, since the different integrals coincide for holomorphic func- 
tions, we shall use indifferently the notation 5 or d. 

As a consequence of Ito's formula, we can deduce the following. 

Proposition 7. Let f : C — > C be a holomorphic function, then 

E ( (im J* f(B s ) dB.} ) = e( (Re jT* f(B s ) dB s 

and 

E^m J* f(B s )dB s } (Re J* f(B s )dB s Y\ =0. 
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Proof. In fact, this is a consequence of Ito's formula. Indeed, let us 
consider a holomorphic function F such that F = f. From Ito's formula, we 
have to show that 

E((Im(F(B t ) - F(z ))) 2 ) = E((Re(F(£ t ) - i^o))) 2 ) 

and 

E((lm(F(B t ) - F(z )))(Re(F(B t ) - F(z )))) = 0. 

Now, notice that after a suitable scaling, these equalities only involve the 
centered reduced normal law, and are, hence, easily checked. □ 

The following proposition provides the skew-product decomposition for 
the fBm. 

Proposition 8 ("Skew-product decomposition for the fBm"). We have 
( t* 5B S \ f f l dB s \ 

Bt=2oexp U^h* oexp U^> 



Proof. If we apply Ito's formula to 

rt dBs 
B s 

we see that it is constant, and so equal to Zq. □ 



Rem / ^\\=ni - :■„• - i-o, 



Notice the important fact that the divergence integral J is local. 
Indeed, if not, we would have 

rt SB S \\ m ( ft B\5B] + B*8B* 
and from Ito's formula, 

E(ln|St|)=E(hi|«b|), 

which is clearly absurd. 

3. Polar decomposition. First of all, we study the polarity of points for 
the complex fractional Brownian motion B issued from zq ^ 0. For this, it 
is enough to adapt the original proof by [6] of the polarity of points for the 
planar Brownian motion. 

Proposition 9. If H > I, then the points are polar for B, that is, 

¥(3t>0,B t = x) = 0. 
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Proof. We proceed in several steps. We assume H > ^ because the 
case H = ^, which corresponds to the case of the Brownian motion, is well 
known. 

In [19] it has proved that a.s. the Hausdorff measure /i^ of the fractional 
Brownian curve is finite, where f{e) = e l l H loglog(l/e). This clearly implies 
that the Lebesgue measure of the fractional Brownian curve is zero. For the 
sake of completeness, we provide a simple and direct proof of this result. Let 
us denote by [i the Lebesgue measure on the plane and note that 

E(n{B t (uj),0<t < 1}) <ttE( sup \B t \ 2 ). 

\o<t<i / 

Hence, 

E(n{B t (u),0<t<l}) <+oo. 
Now, from the scaling property of the fBm, for T £ N*, 

E(fi{B t (uj), <t<T})= T 2H E(fi{B t (uj),0 < t < 1}). 
But, since the increments are stationary, 

T 

E(//{5 t (o;),0 < t < T}) < ^E(^{B t (w),i - 1 < t < i}) 

i=l 

<TE{p{B t (cj),0<t<l}). 

We deduce that 

T 2H E(p{B t (u),0<t < 1}) <TE{^{B t (u),0<t< 1}), 
which implies 

E(fi{B t (io),t>0}) = 0. 
Now, if there exists at least one x £ C \ {^o} such that 

¥(3t>0,B t = x) >0, 
then we claim that, for all z G C \ {^o}, 

F(3t >0,B t = z) >0. 

This follows from the rotational invariance and the scaling property of the 
fBm. Since 

E(p{Bt(u),t> 0})=0, 
we have, for all z S C \ {zo}, 

¥(3t >0,B t = z) = 0. 
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To conclude the proof, we have now to show that 

F(3t>0,B t = zq) =0. 

For this, notice that for, h > 0, the process {B t+ h — Bh)t>o 1S still a frac- 
tional Brownian motion (this is a Gaussian process with the same covariance 
function as B). Hence, if we had 

¥(3t>0,B t = z ) >0, 

then, we would have, for h > 0, 

F{3t>0,B t+h -B h = z )>0, 

which is absurd because zq ^ 0. □ 

Remark 10. It is known (see, e.g., [20]) that, for H G (0, \), the sam- 
ple paths of B have interior points, so that K(p{Bt(u>),0 <t< 1}) > 0. In 
particular, the points are not polar. Notice also that, for ffefOi], B is 
recurrent, whereas it is transient for H > i. 

Assume now that zq = 0. Since is polar for the complex fBm, we can 
consider for t > the radial part 

Pt = \B t \ 

and the angular part 

e, * 



\B, 



Proposition 11. For each t>0,Q t is independent of the process (p s ) s >o 
and uniformly distributed on the unit circle. 

Proof. Let X be an 7£oo-measurable random variable, where 1Z is the 
natural filtration of (p s ) s >o an d let / be a Borel function on the circle. We 
know that the law of B is rotational invariant, and this implies that, for all 
6£ [0,2tt], 

E(Xf(e t ))=E(Xf(e ie 9 t )). 
Integrating this relation yields 

E(X/(G t )) = E (x± fj f(e t9 Q t ) de^j . 

Since the Lebesgue measure on the circle is invariant by translation, we 
deduce that 

1 /" 27r 
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is constant and does not depend on Qf This implies that 

E(xf(e t ))=E(x)±- r f(e i9 )d9, 

2ir Jo 

which gives the expected result. □ 

4. Asymptotics of the skew-product and limit theorems. We come back 
to the case zq ^ and study the asymptotics of J ^p 4 - = f as t tends 
to infinity. First we present some preliminary results. 

Lemma 12. For any real number k > 1, we have 

<- k alB, 



E 



i B s - z 



< +oo. 



Proof. In order to show this property, we express the path-wise Riemann- 
Stieltjes integral in terms of fractional derivatives, using the fractional inte- 
gration by parts formula (see [22]): 

r/R r k / 1 \ 

D i+ (-5— T )(s)Dlz a B k _(s)ds, 



i B s -z Ji + \B.-zo 

where 1 — H < a < h , and the left and right fractional derivatives are given 
by 



ni-«n f \ -1 ( B k -B s f k B y - B s 

Di Bk-(s) = —, , — rz ha / t — - — tts — ay 

and 



{(Bs-zor 1 +a r VCgg - go) - i/(fl y - go) ^ 



r(i-a)V (s-i) a A (s-y) aH 

Clearly, for all p > 1 , 



(4.1) E sup |D^ a S fc _( fl )|f <oo. 

\Ks<k 



On the other hand, for 1 < p < 2 

'i ^B7-z | p (s-l) a?; 
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and 



1 VI 



(4.3) 



/" 7 \ l By Bs \ t r dy] P ds 

I J x ( s -y)*+i\B s -zo\\By-z \ u ) 



< k p E 



pk rs ^ 

G JiJi (s - y)P^+ 1 ~H+e)\B s - z Q \P\By - z \p dydS . 



(G™')) 1/q ' sup (e( ,„ 1 — ^ 29 <oo, 



where 



i< s <fcV \\B s -z \ 2 p<i 

By — B. 



G = sup 



i<y<s<k (s - y) H £ 

g + q> 

(4.2) and (4.3) imply the desired result. □ 



provided | + ^ = 1, 2pq < 2, and p(a + 1 - H + e) < 1. The estimates (4.1), 



Remark 13. The previous result is not true for H = |, that is, in the 
case of the planar Brownian motion. 

Now, we prove the ergodicity of the scaling for the fractional Brownian 
motion. 

Lemma 14. For k £ \ {1}, the transformation on the path space 

Tk : uj — > ^pr which preserves the law of a one- dimensional fBm issued from 
is ergodic. 

Proof. Let (3 be a one-dimensional fractional Brownian motion issued 
from 0. It suffices to show that is mixing. For this, we note that, for 
s,t > and n G N*, 



E (o-fS)^ 



{yS ^ + k 2nH t 2H_ ls _ k n tl 2 H] 



so that 



n— »+oo \ k 

which implies the mixing property of T^. □ 

Before we turn to the skew-product decomposition, we state a last pre- 
liminary result on the -^-variation of the integral Jq We refer to [5] for 
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similar results on the -^-variation of divergence integrals with application 
to fractional Bessel processes. 

Proposition 15. Set U = for i = 0,1, ... ,n. Then 



n-1 

E 

i=0 



dB< 



l/H 



ds 



almost surely, as n tends to infinity, where c = E{\B\\ l / H ) 



Proof. For any complex- valued process ut whose paths are locally 7- 
Holder continuous, with 7 > 1 — H, we set 



n— 1 »+. . 



l/H 



i=0 

Then, by the Holder inequality, we obtain 

\V n (u)-V n (v)\ <tfV( / 1 (u s -v s )dB 
H - =0 V Ju 



(4.4) 



ti+l 



u s dB s 



l/H 


fU+l 




+ 


/ u s <iB s 






Jti 





< -\V n (u - v)\ H {V n {uf- M + V n {v) l - H ). 



If u t is a step process, then, from the properties of the fBm (see [15]), we 
know that V n (u) converges almost surely to c Jq {ugl 1 ^ 11 ds. We can write 



Vn 



ds 



\B 8 \V H 



< 



Vn 



Vn(Um) 

u m (s)\ l ' H ds 

Q"n,m ~\~ b n ,m ~\~ Cjni 



+ 



V n (u m )-c \u m (s)\ 1/H ds 



ds 



\B x \y H 



where 



m— 1 



1 



Um(t)= J2 R 1 (t j - 1 ,t j ]^- 

3=0 ^J-i 

Using (4.4), we can estimate the first summand as follows: 



Q"n,m — 



II 



v„ 



B 



H 



Vn 



B 



1-H 



+ V n (u r 



A-H 
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1 

< — 
- H 



Vax l/H [ ^ ( — - u m ] dB s 



H 



x(v&r 1/H ^ ~Q dB ^j +Vari /^(^ u mdB, 

where Var 1 /^ denotes the total variation of order jj on the interval [0,t]. 
Then Love- Young's inequality (see [21]) implies that lim m sup n a n ^ m = 0. 
From the previous results, we know that, for any fixed m, lim n b n ^ m = 0, 
and clearly, lim m c m = 0. This completes the proof. □ 

We can now turn to the asymptotics of the skew-product decomposition. 
The main ingredient of what follows is the ergodic theorem. 



Proposition 16. We have 



1 /•* SB 



logt JO B s i^+oo 

and 

1 I* ds a . s . r(l-l/(2il)) 



logt Jo \B s \ l l H t-,+oo 2V(2H) 
Proof. We first write 

B t = zo + Pt, 

where (3 is an fBm started at 0. 
The key point is that 

d(3 s ds 
and 



are invariant by scaling and this suggests to use Birkhoff-Khintchine's er- 
godic theorem. 

Let k > 0, since the law of the standard fBm is invariant under the trans- 
formation Tk : to — > , by Birkhoff-Khintchine's theorem and Lemmas 12 
and 14, we have 

1 ^ rk dP k n s /k nH a . s . 

where Y is constant. Hence, 
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We conclude that 

l r kN dp, 

N J i (3 S N^+oo 

Since k was arbitrary, 



1 /■* dp. 



log til (3s t^+oo 

Now, because B is transient, 

1 fdPs 1 [ l dB s 



log* J i (3 S logt Jo -B s ' 

This implies 

1 



>s a.s. y 



logWo B s t^+oo 

Let us now show that 

Re(Y) = H. 

For this, we must check that 

e(«) ^ if, 

V log t / t^+oo 

which is a direct consequence of the scaling property for (5. We deduce that 

ReY = H. 

Now, we note that 

ImY = 

because of the rotational invariance of the fBm. 

Since the second limit is obtained exactly by the same way, we omit the 
details of its proof. Let us just mention that it is based on the following 
computation: 



t->+oo VlogUo IBs] 1 / 11 J VliVl 1 /^ 

where N denotes a random variable with standard normal distribution on 
E 2 . □ 

Remark 17. Unlike what happens in the case of the standard Brownian 
motion, the process (J ^^-)t>o is not a complex fBm P time-changed with 
the clock 

ds 



o 
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Indeed, otherwise, we would have 

_L_ f t S_B 1= Pj*ds/\B a \^ H Q 
log t Jo B s log t t^+oo 

because 

0t a,s,^ q 

t t—* + OD 

Remark 18. It is also possible to apply the ergodic theorem to obtain, 
more generally, jointly with Proposition 11, the following limit theorem: 

1 f« 1 JP.\, a.s. r(i-i/(2g)) f 
log* A W/ hJ \\Ps\) ^ 2V(2H) J S J^ 

where / is an integrable function on the unit circle § endowed with its 
normalized Haar measure a and j3 is a two-dimensional fBm issued from 0. 

As a direct corollary of this, and by time inversion we, deduce the following 
corollary. 

Corollary 19. Let (Bt)t>o be a complex fBm with Hurst parameter 
H > i and started at 0. We have, for any z G C, 



logl-Bt + d a .s. TT , log\B t + zt 



2H\ 



H and — - ' H. 

logt t-^+cxD logt t^0+ 

Proof. The first part of the proposition stems immediately from the 
skew-product decomposition, indeed, 

log |4| 1 ^(f'SB s 



log* logt \Jo B s 
For the second part, we use a time inversion argument. Indeed, we can write 

B t = t 2H p 1 / t , 

where f3 is a planar fractional Brownian motion issued from 0, and then we 
apply the first part of the proposition to /3. □ 

Observe that, from the law of iterated logarithm (see, e.g., [9] and [10]), 
it is known that, for any H £ (0, 1), the following inequality holds: 



\B t \ <A(l + t H )^/lnln(3 + t). 
Before we conclude this section, we give a small extension of Proposition 16. 
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Proposition 20. Let (Bt)t>o be a complex fBm with Hurst parameter 
H > 2 and started at zq j^O. Let f : C — ► C be a holomorphic function on the 
whole plane such that /(0) = 0. Then 



L [ t f(±) 5Bs *±> H f'(0). 
>gtj \B S J s t-+oo v ; 



log 

Proof. Indeed, we can write for z £ C* and some holomorphic function 

V Z J z z z \z 



It stems from this that 

logtJ J \B s J" 6 logt Jo B s ' logtj B 2 * \B, 
Now, thanks to Ito's formula, 



1 '* VteW- 1 



i \ / i 



9{ V )-\B t 



logtJo Bf \B S ) log* 
which allows to conclude according to Proposition 16. □ 

5. Study of the windings. In this section we shall again assume that 
Zq 0. We denote 

1 B 2 S dBl - Bl dB 2 

From the skew-product decomposition of the two-dimensional fractional 
Brownian motion, the process 9 characterizes the windings of B. Though 
the complete study of the asymptotics of 9 is not yet well understood, we 
present some results. First, it is easy to study the windings of B up to 2n. 



Proposition 21. We have 



1 



Me inet ) ~ 

V ' n->+oo t nH 

Proof. From the skew-product decomposition, we get 

T ± = exp(i9 t ), t>0. 
\Bt\ 

Hence, for n £ N, 

■(t H N + z y 



E(exp(m0 t )) =E 



\t H N + z \ n J' 
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where N is a two-dimensional standard normal law. The result follows then 
from a straightforward computation on the two-dimensional standard nor- 
mal law. □ 



We conclude now the paper with a limit theorem for the functional of the 
two-dimensional fractional Brownian motion 



(5.1) 



* B*dB\-B\dBl 

S 2H 



which looks like the windings. Formula (2.3) allows us to compute the vari- 
ance of this process: 

E(Z?)=2a H / (rs)~ 2H R(s,r)\r- s\ 2H ~ 2 drds 



l Jl 



(5.2) -2a 2 H I (sr)- 2H l, e<ra<s} \r-a\ 2H - 2 \9-s\ 2H - 2 drdsdeda 
J[i,t] 4 ~ ' ~ 

<C H logt. 



Proposition 22. Let Z t be defined by (5.1). Then -j=k = converges in 
law, as t tends to infinity, to the normal distribution N(0,a 2 ), where 

a 2 = 4H(2H-1)£ (y~ 2H R(y, 1) - H log y' 1 - Hp (3^-) ) (1 - v) 2H ~ 2 dy 
and 

fi 

\2H-l/ . . \-2H 



0(y)= / (l- X y^{x + y)-^dx. 
Jo 

Proof. Define, for A G R and t>l, 

At = exp i\- 



and set <p(X, t) = E(A^ ) . Using the change of variable formula for the Riemann- 
Stieltjes integral we can write, for 1 < to < t, 

X t = X tQ + f -^=X u u- 2H (B 2 u dBl - B l u dB 2 u ) 



(5.3) 



Vlogu 

X u Z u du. 



t 2u(logii) 3 / 2 

Using (2.4), we can transform this pathwise integral into a divergence inte- 
gral, plus a complementary term: 

f ~^=X u u~ 2H (Bl dB l u - Bi dB 2 " 

Jtn Vlo, 



ylogu 
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(5.4) = I* -i±=X u u-™(Bl 5Bl - Bl 5B 2 U ) 

Jtn VlogU 



rt n~2H f"U 

+ a H i\ / - / {BlDlX u - BlD 2 r X u )(u - r) 2H ~ 2 drdu. 
Jtn V loe u Jo 



ydog u Jo 

Substituting (5.4) into (5.3) and taking expectations, we obtain 
(p(X,t) = <p(\,t ) 

-a H \ 2 /*t E r X u (B 2 DlZ u -BlD 2 Z u )(u-r) 2H ~ 2 drdu 

J tn logn Jo 



to logw 
* A dip 
to 



(A, u) du. 



2u log u d\ 

Now we apply the duality relationship (2.2) and we get 
E(X u BlD l r Z u ) = K(X u {D 2 DlZ Ul l (0 , tt )) w ) 

= ^{X u )(D 2 DlZ u , 1(o,«))h 
iX 



+ 



-=E(X M ( J D 2 Z U , 1 (0)U) ) W £>*Z U ) 
ylogu 



and we obtain 

¥>(A,i) = 9?(A,i ) 



/•t n~2H 

«^A 2 / E(X U 

Jfa, logU 



'to 

X 



logu 

/ ((D 2 DlZ u ,l( ^) n 

J 

ln2 7 i \ \t.. „\2H-2 



- {D^D z r Z u , 1 {0 ,u))h)( u ~ rY H - z drdu 
Jtn 2u log U OA 



where 

i/)(\,t) = —i\ 3 aH 



t U -2H 



to (logu) 3 / 2 

- (D}Z U , l {0 , u) ) H D 2 Z u )](u - r) 2H ~ 2 drdu. 
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We have, for r <u, 

DlZ u 



b 2 r 


SB 2 




Q2H 


f u SB] 




Jr Q 2H 





Hence, 

(D}D 2 Z U , l(o,it))w = (^~ 2Hl (r,n) - r ~ 2Hl (0,r)> 1(0,m)} 



H 



u ru 

«- 2 ^a-6\ 2H - 2 d6do -r- 2H R(r,> 



Jr 



H [ U e~ 2H (6 2H ~ l + {u- Of 11 " 1 ) d6 - r~ 2H R(r, u) 



11 

H\og- + H(3 

r 



(—)-r- 2H R(r,u), 
\u — r J 



where 



and, similarly, 



0(y)= f\l-x) 2H - l (x + y y 2H dx, 
Jo 



{D 2 DlZ u , 1(o,«)) w — ( r 2i? l(o,r) ~ O 2H l(r,u), 1(o,u))h 

= r- 2H R(r, u) - fflog - - Hp( — — 

r \u — r 

Hence, 



ip(\,t) = v(\,t ) 

- 2a H X 2 [ y E(X U ) 

J to logn 



'*o 

X 



logu 

f ( r' 2H R(r, u)-Hlog--Hp( —) 
Jo \ r \u — rj 

(5.5) x {u-r) 2H - 2 drdu 

rt X dy 



to 



2u log u dX 



(A, u) du + ip(X, t) 



f 1 

<p(X,t ) - 2a H X 2 / — : <p(\,u)p(u)du 

r* A dip 
1 1 2ulogu d\ 



f f X dp 
J to 2u log u dX 
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where 
P{z) 



f xz [y~ 2H R{y, i) - i^iogy" 1 - ^)) (l - y) 2H ~ 2 dy. 



We claim that lim^oo ip(X, t) = 0. In order to show this property, let us write 
E(X U (D 2 Z U , l^ u )) H D r Z u ) 

Bl rdBl 



E ^X u (D 2 Z u , 1 ( 0iU ) ) H f - 



(DfX u ,9 2 l( r , u ))w)) 



iA 



E(X„(D^ tt) l (0 , u )) H (r- 2H ( J D 2 Z it , l (0jr) ) w 



{DfZ u ,9 2H l^ u )) n )). 



Therefore, 

•0(A, i) = — i\ 3 ctH 
x E 



* u~ 2H 



:X n , 



to 



log m 3 / 2 u 

ru 

x / ((D 2 Z u ,l(o,«))w 

<J 

x (r~ 2H (D 2 Z U , 1(0,1-)) « - (D 2 Z u ,9~ 2H l^ r ^) n j 
+ (D}Z U , l(o,n))« 

x (r- 2 ^( J D!Z u ,l (0ir) ) w - (D. 1 ^^- 2if l(r,n))K)) 
x (u — r) 2H ~ 2 drdu. 

We have 

(D}Z U , 1( ,u))h = olh 



o io V 



Hence, 

IKAZi, l(0,u))wll2 - a # 



JO 



-2H 



Bi 



r-9\ ltl ~ l drd9. 



5Bl 



Q2H 



\r-9\ 2H ~ 2 drd9 



<C r~ H \r-9\ ZM - 2 drd9 
Jo Jo 

<Cu H , 
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because, applying the inequality proved in [7], we obtain, for some constant 

CH, 

[ u 5B% H/3-2H-1 II <- „ „-H 

-Q2H ^<C H \W 1 (r,«)llLl/ff(R + )< C H r . 

Similarly, 

(D}Z u ,9~ 2H l( ru )) H = an 



u f u ( B 



r W H J a 9 2H J 



SBl 



\a-0 



2H-2a-2H 



d6da 



and 



(D}Z u ,0- 2H l {r:U) ) n \\ 2 < CJ J a-"\a-e\ z "-' z 9-' ztl dOda 
= Cdr~ H , 

d = J J y~ H \y-x\ 2H - 2 x- 2H dxdy <oo. 



-Hi q\2H-2q-2H 



where 



Also, 



(D}Z U , l(o, r ))wll2 - a H 



o Jo 

u rr 



5B 2 Bl 



a 9 2H a 2H 



\a-e\ 2H - 2 d6da 



<C a- H \a-9\ 2H - 2 d9da 
Jo Jo 

<Cr H . 



Finally, 



(5.6) 



m\,t)\<c 



< 



c 



where the constant C depends on to, H and A. Similarly, we can derive the 
estimate 



(5.7) 

Equation (5.5) yields 

dip 



< 



C 



t{\ogtfl 2 ' 



.-2a H X 2 ^ A_9^ 9^ 

dt tlogt 2t\ogtd\ dt 
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Hence, 

J to I Jto ulogu JV 2s logs OA ds J 

The integral exp{— f* Q —j^j-^-ds} behaves as j^p- as t tends to infinity. 
We have, from (5.2), 

dtp 



(5.9) 



dX 



< -P=\E(Z t X t )\<\\\^. 
Vlogi 



Then, (5.9) and (5.7) imply that lim^oo ip(\,t) = <p(\,oo) exists. By a sim- 
ilar argument, differentiating equation (5.8) with respect to A, we can show 
that limt—nx, ^ = ^(A,oo) exists. We claim that 

(5.10) -2a H Ap(oo)99(A,oo) - ^ (A, oo) = 0. 

Otherwise, ^ behaves as t ^ t when t tends to infinity, where C is a constant 
different from zero, which is contradictory because we would have 

f(X,t) - Clog log t. 

t—>+oo 

Equality (5.10) implies that 

and Si converges in distribution as t tends to infinity to the normal 
distribution N(0, 4o#p(oo)). □ 

Notice also that, as H tends to h, the variance a 2 in the preceding propo- 
sition converges to 2. This is reasonable, because if Bt is a standard two- 
dimensional Brownian motion, , L —* — s 2 — *- converges in distribu- 

tion as t tends to infinity to the normal law iV(0,2). Indeed, the Girsanov 
theorem implies that 

/ / iX [* B 2 S dBl- Bl dB 2 s X 2 ft \B S \ 2 , \\ 

1 rt \B 2 

and, applying the ergodic theorem, we have that 2 \ogt fi J ^ s conver g es 
almost surely to one. Nevertheless, we could expect that, in the case of the 
fBm which is transient, the functional 

* B 2 dB]- B]dB 2 s 



S 2H 
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is closer to the windings of B than it is the case for the planar Brownian 
motion for which it is well known, since [18], that they behave like ^Clogi, 
where C is a Cauchy law with parameter 1. 
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